Quasi Non-linear Evolution of Stochastic Bias 



Atsushi Taruya 1 , Kazuya Koyama 2 and Jiro Soda 1 
1, Department of Fundamental Sciences, FIHS, Kyoto University, Kyoto, 606-8501, Japan 
2, Graduate School of Human and Environment Studies, Kyoto University, Kyoto, 606-8501, 

Japan 

ABSTRACT 

It is generally believed that the spatial distribution of galaxies does not trace that of 
the total mass. The understanding of the bias effect is therefore necessary to determine 
the cosmological parameters and the primordial density fluctuation spectrum from the 
galaxy survey. The deterministic description of bias may not be appropriate because 
of the various stochasticity of galaxy formation process. In nature, the biasing is 
epoch dependent and recent deep survey of the galaxy shows the large biasing at high 
redshift. Hence, we investigate quasi non-linear evolution of the stochastic bias by 
using the tree level perturbation method. Especially, the influence of the initial cross 
correlation on the evolution of the skewness and the bi-spectrum is examined in detail. 
We find that the non-linear bias can be generated dynamically. The small value of 
the initial cross correlation can bend the 5 g -5 m relation effectively and easily lead to 
the negative curvature (62 < 0). We also propose a method to predict the bias, cross 
correlation and skewness at high redshift. As an illustration, the possibility of the 
large biasing at high redshift is discussed. Provided the present bias parameter as 
b = 1.5 and Q, = 1.0, we predict the large scale bias as b = 4.63 at z = 3 by fitting the 
bi-spectrum to the Lick catalog data. Our results will be important for the future deep 
sky survey. 

Subject headings: cosmology:theory— large scale structure of universe 



1. Introduction 



The large scale structure in the universe has evolved from the primordial mass fluctuations 
according to the gravitational instability. In standard picture of structure formation, the density 
fluctuations are produced by the quantum fluctuation during the inflationary stage of the early 
universe. The CMB anisotropy observed by COBE reflects such fluctuations at the last scattering 
surface. Subsequently, the growth of the mass fluctuation leads to formation of galaxies and 
evolution of galaxy distribution which is observed by the galaxy survey such as APM, CfA survey, 
etc. The most important issue in modern cosmology is to construct a consistent history of our 
universe from inflation to galaxy clustering. 
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Statistical quantities of large scale structure provide us many important information to 
understand the evolution of our universe. From the power spectrum or the two-point correlation 
function of galaxies, we can evaluate the density parameter fig by comparing the observations in 
the redshift space with those in the real space ( [Hamilton 1997 ). In general, the n-point galaxy 



correlation functions characterize the large scale structure of the universe. These quantities play a 
significant role to test the prediction of inflation theory. However, there exists an uncertainty in 
the relation between the distribution of galaxy and that of total mass in the universe. Galaxies 
have been formed at the high density regions of the mass fluctuations. Due to the lack of our 
knowledge about galaxy formation process, we do not have the definite answer how much fraction 
of the total mass is the luminous matter. This uncertainty is quantified by the bias parameter of 
the galaxy distribution, which significantly affects the determination of cosmological parameters 
and the tests of inflation theory. 

To proceed further, it is necessary to model the galaxy distribution in the given mass 
distribution. The linear bias gives a simple relation between the galaxy and mass when both 
fluctuations are small enough. Denoting the fluctuation of total mass as 8^ and that of the galaxy 
mass as 5 g , we have 

Sg = b S m , (1) 

where b is the linear bias parameter. The assumption ([!]) yields a degeneracy of the parameters, 
(3(Qq) = f2[j' 64 /&, f° r determining Qq ( [Hamilton 1997|) . We need another observation in order to 
lift up this degeneracy. Moreover, we should keep in mind that the relation (|l]) cannot be static. 
Indeed, the bias parameter changes in time by gravitational force even if the galaxy formation 
turned off. The recent observation of the high redshift galaxy survey suggests that the galaxies at 
z ~ 3 are largely biased, whose bias parameter is evaluated as b ~ 6, even if we do not require the 
bias at present ( peacock 1998 ). It is necessary to consider the bias effect by taking into account 
the time evolution. 

The linear bias is appropriate as long as the linear theory is a good approximation to the 
observation. A naive extension of biasing to the non-linear regime is the deterministic non-linear 
bias. Suppose that the galaxy distribution 5 g is a local function of 5 m , we express 8 g in powers of 
8 m ( [Fry fc Gaztahaga 1993|) : 

s s = /(<W = £^(<W n (2) 

n 

The non-linear bias parameter b n appears when we compare the higher order statistics of the 
galaxy distribution with that of the mass fluctuation such as the skewness (S3) and the kurtosis 
(S 4 ) defined by S 3 = (S 3 )/(5 2 ) 2 and S A = ({5 4 ) - 3{5 2 ) 2 )/ {5 2 ) 2 , respectively. The relation becomes 



S ^ = 7£x =b- 1 (S 3 , m + 3^) + 0((5 2 m )). (3) 



S., 9 = ^ {6 2 } f 9)2 = b- 2 (S 4 , m + 12 b fS 3 , m + 1 2 | + 4) + 0((S 2 J), (4) 
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where we put b\ = b. Subscript g and m denote the statistics of the galaxies and the mass 
fluctuation, respectively. Other than the above one-point functions, we can calculate the amplitude 
of the galaxy three point correlation function (bi-spectrum) Q g which also has the relation 
Qg = Qm/b + b2/b 2 ( [Fry 1994 ). Since the statistical quantities for mass distribution can be 



obtained analytically under the assumption of gravitational instability, we can determine the bias 
parameters b, 62, 63, • • • , by combining the observation of S 9t 3, S g ^ and with that of Q g flFriemari 



fc Gaztahaga 1994 , Gaztahaga fc Frieman 1994| ). As for the time evolution of the bias parameter, 



Fry has discussed its influence on the skewness and the bi-spectrum ( Fry 1996 ). He found that 



5*3^ and Q g asymptotically approach the constant values of S3 )Tn and Q m , respectively. 

There has been some confusions for the galactic bias effect since Kaiser provided a simple 
bias mechanism ( Kaiser T98|). His original idea is to explain the enhancement of the two 
point correlation function of the rich clusters using the statistical feature of galaxy distribution. 
Assuming the rich clusters are formed where the fluctuation 5 g averaged over a larger scale exceeds 
a certain threshold, he derived a simple relation between the two-point correlation function of 
galaxies and that of rich clusters. It should be stressed that the bias parameter introduced by 
Kaiser is not defined by density fluctuations but correlation functions. In our case, the galaxy 
correlation function may relate with the mass correlation function, however, the deterministic 
relation between the variables 5 g and 5 m does not necessarily follow from the relation between 
correlation functions. 

According to the inflationary scenario for the structure formation, the primordial density 
fluctuation is generated quantum mechanically and it is likely to have stochastic nature. ( [Linde 



1990). Recently, the galaxy formation process is recognized as non-linear and stochastic one (Cen 



& Ostriker 1992). Therefore we should treat both S g and S m as independent stochastic variables. 



The stochastic treatment of the bias effect is referred to as the stochastic bias, which has been 



introduced by Dekel and Lahav (Dekel 1997, Lahav 1996). In stochastic bias, the relation between 



the galaxy and the total mass can be described by the correlation functions. When the fluctuations 
are small enough and the linear perturbation is valid, we need three parameters defined below: 

a 2 = (5 2 ) b 2 = ^- r - {5m5g) (5) 

The important ingredient of the stochastic bias is the cross correlation r which is absent in the 
simple relation of linear bias ([!]). This is an extension of the Kaiser's bias prescription. Matarrese 
et al. shows that the cross correlation between the regions with different thresholds of density 
peak naturally arises in general situation including both Gaussian and non-Gaussian spatial 
distribution ( Matarrese et a/.1986j ). Mo & White discussed the spatial distribution of galactic 



halloes using the Press & Schechter formalism and evaluated the cross correlation by comparison 
with numerical simulation ( [Mo Sz White 19*96]) . Recently, Catelan et al. extended this formalism 



(Catelan et al. 1998a, Catelan et al. 1998b, Porciani et al. 1998). We think that the cross correlation 



plays an important role to describe the galactic bias effect. 
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A naive introduction of cross correlation increases the unknown parameters and may 
complicate the determination of the cosmological parameters from the measurement of galaxy 
survey. Pen has discussed how to determine the parameters b, r, a and the other non-Gaussian 
variables. He concluded that the observation of redshift space distortion is useful to understand 
these parameters and the density parameter Slo ( pen 1997| ). However he did not consider the time 
evolution. On the other hand, Tegmark & Peebles studied the linear evolution of the stochastic 



bias after the galaxy formation epoch and on-going the galaxy formation ( [Tegmark fc Peebles 



1998 ). They found that the bias parameter b and the cross correlation r approaches unity even if 



galaxy formation never ends. However, their analysis is restricted to the linear evolution. 

Purpose of the present paper is to understand the quasi non-linear time evolution of 
the stochastic bias by observing the three-point correlation functions ( [Fry 1994| ). The future 
development of galaxy survey such as Sloan Digital Sky Survey(SDSS) and two degree field survey 
(2dF) will provide us an enormous data of the large scale structure. In the future, the time 
evolution of galaxy distribution will be observed. Study of the bias evolution is of importance to 
separate the effect of the time evolution due to the galaxy formation process and the growth due 
to gravity. To investigate this, we must understand the qualitative and/or quantitative behavior 
of higher order statistics under the influence of gravity and recognize a role of the initial cross 
correlation which is important quantity in the stochastic bias. 

In this paper, we shall pay an attention to the time evolution of skewness and bi-spectrum 
in the stochastic description of the bias effect. We discuss the basic formalism of stochastic bias 
and derive the evolution equations for galaxy and mass distribution in Sec.||. Based on this 
formulation, we first evaluate the bias parameter and the cross correlation at the tree level. In 
Sec.||, we develop the second order perturbation and study the time evolution of three-point 
function. The effect of the initial cross correlation is investigated by comparing with the various 
choices of the parameters. We find that the initial cross correlation significantly affects the 
non-linear relation between 5 m and 5 g which differs from the deterministic bias. Sec|I| is devoted 
to the discussion about the observation from the point of the stochastic biasing. Using the Lick 
catalog data and the present bias parameter inferred by the IRAS galaxy survey, we shall show 
that the high redshift galaxies are strongly biased, which will be observed by the future galaxy 
survey. We summarize the results briefly and present the future prospect in Sec.|5[ 



2. Basic Formulation of Stochastic Bias 

This section is devoted to a basic formulation of stochastic bias for further analysis of next 
section. We first consider the stochastic description of fluctuation for the galaxy distribution and 
the total mass distribution. In the case of weakly non-linear evolution, the generating functional 
for the stochastic variables of galaxy and total mass is constructed. Subsequently, we give the 
evolution equations and initial conditions. We will analyze these equations perturbatively and 
discuss the linear evolution of the stochastic bias. 
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2.1. Stochastic description 



The basic quantities to treat the stochastic bias are the fluctuations of total mass and the 
galaxy distribution, 5 m and S g . They are defined by 



S m (x) 



p{x) ~ P 

— 1 

P 



5g(x) 



n g (x) - n g 



(6) 



where barred quantities means homogeneous averaged density. We shall investigate the skewness 
and the bi-spectrum of the galaxy distribution 5 g . We evaluate them using the smoothed density 
field by the window function: 



6m, g (R) = d 3 x- W R {x)5 myg (x). 



(7) 



The window function Wr{x) we adopt here is the spherical top-hat smoothing with the radius R 
defined bellow: 



W R {x) 



1 (|x| < R), 
( x > R). 



(8) 



As we are interested in the time evolution of bias effect under the influence of gravity, we 
ignore the galaxy formation process in consideration of the evolution of galaxy distribution. The 
stochastic bias is described by the stochastic variables 5 m and S g . Their stochasticity is determined 
by the probability distribution functional V\8 m {x), S g (x)]. Once V[5 m (x), 5 g (x)] is given, all the 
information of the distribution functional is encoded in the correlation functions of 5 m and S g , and 
all of the correlation functions can be obtained from the following generating functional: 

Z[J m ,J g ] = J J V5 m {x)V5 g (x)-V[5 

m (x) ) \X)\ 

d 3 xW R {x) {J m 5 m (x) + J g 5 g (x)} 



x exp 

where J m , J g are external source terms. Using (|9|), the n-th moment of the filtered one-point 
function is written by 



(9) 



({S m (R)} j ■ {5 9 (R)} n - j ) = i~ n —] 



5 n Z 



(10) 



Jm — J a — 



The probability functional V[5 m (x),S g (x)] can be determined in principle by the study of 
galaxy formation and evolution of fluctuations before galaxy formation. In our prescription of 
stochastic bias, we give V[5 m (x),S g (x)] as a parameterized function. A simple choice is Gaussian: 



V[S m (x),5 g (x)] = AT 1 exp 



-(8 m ,8 g )G 1 



(11) 



- 6 - 



where G is a 2 x 2 matrix and TV is a normalization constant. The non- vanishing off-diagonal 
component of G represents the cross correlation between S m and 5 g . We obtain 



(S m (R)S g (R)) = / d 6 xd 6 yW R {x)W R {y){G} 12 5 m {x)5 g (x) 



(12) 



To take into account the non-linearity of gravitational evolution and the formation process of 
galaxies, the Gaussian statistics does not lead to a correct description. Pen has discussed the 



deviation from the Gaussian distribution by using Edgeworth expansion( Juszkiewicz et al. 1995 ) 



although his analysis is restricted to the static case ( Pen 1997 ) 



Since the density fluctuations are small on large scales, perturbative analysis is applicable to 
study the evolution of density fluctuation. In this case, the non-Gaussian stochastic distribution 
for 5 m and 5 g can be expressed in powers of the Gaussian variables A m and A g whose variances 
are small enough. We shall write the initial condition at the end of galaxy formation (t = ij) as 



/(Am), 



2(A 9 



(13) 



If we know the function / and g and stochastic property for A m and A g , the initial condition for 
5 m and 5 g is determined. The time evolution may change the initial statistics. We can obtain the 
correlation functions after the galaxy formation epoch by solving the evolution equations for 5 m 
and 6 g . Then the time-dependent generating functional becomes 



-2^ [Jm , Jg ! 1\ 



VA m {x)VA g (x) ■ N' 1 exp 



-(A m , A g )G 



-l ( A„ 
A c 



x exp 



i J d 3 xW R (x){J m 5 m (A m ,A g ,t) + J g 5g(A m ,A g ,t)} 



The time dependence of 5 m and 5 g are obtained from the perturbative expansion: 
6 m (A m , A g , t) = 6® + 4 2) + ■ • • , S g (A m , Ag, t) = 4 1 ) + 5^ + ■■■ 



(14) 



(15) 



Since the n-th order perturbed variables 8m]g have the same order of magnitude as (A m>9 ) n , the 
initial condition given by (13) assigns the solutions (|l5|) order by order. From ( |14| ) and (|l5|), the 
lowest order non-vanishing n-th moment can be evaluated from the results up to the (n — l)-th 
order perturbation. It is composed of "tree" diagram with no loop, in graphical representation 



of the perturbation theory ( Fry 1984 ). The tree- level third order moment for galaxy distribution 
becomes 



(MR)? 



5 3 Z 



3([4 1 )(i?)] 2 4 2 )( J R)). 



(16) 



From (|14j), in principle, we can obtain the joint distribution function of 5 m and 5 g by Laplace 
transformation. 
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2.2. Evolution equations and initial condition 

The evolution of the galaxy and total mass distribution after galaxy formation is determined 
by gravity including the effect of cosmic expansion. The total mass density p(x) is approximated 
by a non-relativistic pressureless fluid. The averaged homogeneous part p is obtained by solving 
the FRW equations: 



H 2 = —p+--- ; (K = 0,±1), (17) 



d\ 2 tt2 8vrG A K 

a J 6 6 a z 

a AttG A 

- = — ( 18 ) 

a 3 3 



where a is the expansion factor, K is the spatial curvature, and A is a cosmological constant. 
Using these variables, the density parameter Q is defined by 

Let us consider the fluctuating part of the mass distribution 5 m . It obeys the equation of 
continuity and the peculiar velocity field v is determined by the Euler equation in the presence 
of gravitational potential. On large scales, the assumption that the velocity field is irrotational 
would be valid. Then, the basic quantities for describing the dynamics are reduced to 5 m and 
9 = V • v/(aH). The evolution equations become ( Peebles 1980| ) 

^ + ^+-v-M = o, (20) 

ot a 

d6 ( O A \ 3 1 



As for the galaxy distribution, 6 g should satisfy equation of continuity as long as the galaxy 
formation is not efficient. After formation epoch, the galaxies move along the velocity field 



determined by the gravitational potential. Therefore we have ( Fry 1996 ) 



^k + H e + -V-(8 g v) = 0. (22) 
ot a 



Eqs.(pOD, (plD, and (22) are our basic equations to determine the evolution of 5 m and S g . 

We next consider the initial condition (^) given at the end of galaxy formation t = t%. For 
the total mass fluctuation, we believe that 6 m is produced during the very early stage of the 
universe. In standard scenario of the inflationary universe, the density fluctuation is generated by 
the quantum fluctuation and may have the Gaussian statistics. We regard such fluctuations as 
A m . After the inflation, the gravitational instability incorporates the deviation from the Gaussian 
fluctuation. Since the galaxy formation does not affect the evolution of 5 m on large scales, the 
total mass fluctuation 5 m at the end of galaxy formation can be determined by the gravitational 
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instability. At that time, the growing mode is dominant. Therefore, we obtain the initial condition 
= /(A TO ) from the evolution equations by dropping the decaying mode. 

On the other hand, the fluctuation of galaxy number density is induced by the galaxy 
formation. To determine the function g(A g ), we need to know halo formation processes. An 
analytic model for the spatial clustering of haloes is discussed by Mo & White and several authors 
( Mo &: White 1996| ). Here, we rather treat g(A g ) as a parameterized function, whose unknown 
parameters are determined by the observation of galaxy survey. Assuming g(A g ) as a local 
function of A g , we have 

g(A g ) = A g + ^(A 2 g -(A 2 g )) + .... (23) 



The remaining task is to specify the stochastic property of A m and A g . Since these are the 
Gaussian variables given by (||), their statistics can be characterized completely by the three 
parameters below: 

*o = (A^), &3 = 7^, ro= - (AgAm) xl/2 , (24) 



which is equivalent to determining the matrix G. Notice that — 1 < ro < 1 can be deduced from 
the Schwarz inequality. We simply assume that the parameters bo, rg are constant, which comes 



from the fact that there is no evidence of the scale-dependent bias on large scales ( Mann et al. 



1997 ). Therefore power spectra for the galaxy two-point correlation function and the galaxy-mass 
cross correlation are simply represented by the power spectrum P(k) of the total mass fluctuation 
given by 

r rl 3 k 

* 2 o(R) = J j^W 2 R (kR)P(k), (25) 
where Wn(kR) is the Fourier transform of the top-hat window function, 

W R (kR) = —^-3 [sin (kR) - kR cos (kR)} . (26) 



2.3. Variance and covariance 

We are in a position to study the evolution of stochastic bias. To begin with, it is necessary 
to study the variance and covariance. To calculate these quantity at the tree level, the linear 
perturbation is sufficient. In next section, we shall develop the second order perturbation and 
analyze the skewness and bi-spectrum. 

Initial conditions and the evolution equations (p0|), (|2l|), and (j2^) yield the linear order 
solutions: 

6&\x,t) = A m (x)D+(t), (27) 
6®(x,t) = A m (x)(D+(t)-l)+A g (x), (28) 
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where the function D + (t) denotes the solution of growing mode by setting D + {ti) = 1, which 
satisfies 

D + + 2HD+ - ^H 2 nD + = 0. (29) 
In Einstein-de Sitter universe (Q = 1), we have D + (t) = a(t)/a(ti). 

Since the stochastic property is already given by (|14|) and (p4|), we can obtain the relation 
between 5 m and 5 g after galaxy formation. Similar to the expression (|24|), we define the 
time-dependent parameters as follows: 

^(t) S <&(t)>, *W = ^> = <WM '» (30) 

{dm{t)} ((si(t))(si(t))) 

Substituting @ and @ into @), we get 



/(Z? + -l)* + 2ftor (S + - !)+*>§ 
K*) = J > ( 31 ) 



These equations coincides with the result of Tegmark & Peebles in the absence of galaxy 
formation ( [Tegmark Peebles 1998| ). The time dependent correlation functions ( |3~l|) show that 
the bias parameter 6(t) and the cross correlation r{t) approach unity asymptotically as the 
fluctuations grow. Note that the deterministic bias relation restricts on the parameter ro = 1, 
which leads to r(t) = 1. The value ro = means that 5 m and 6 g have no correlation initially. 
Subsequent evolution of 5 g and 5 m develops the correlations and ultimately makes the complete 
correlation, i.e, r = 1 in the case of Einstein-de Sitter universe. This indicates that the cross 
correlation significantly affects the evolution of the higher order statistics. The functional form of 
the bias is non-linear, in general. This relation gives the linear bias relation (|l|) if we restrict us to 
the linear region, however, this is not sufficient obviously. 



3. Skewness and Bi-spectrum 

Following the formulation of the stochastic bias in previous section, we investigate the bias 
effect on the skewness and the bi-spectrum. First of all, the solutions of second order perturbation 



are given. Using this result and the analysis in Sec. 2. 3, we evaluate the time evolution of the 



skewness and bi-spectrum of galaxies and analyze the influence of cross correlation r on them in 



Sec. 3.1 and Sec.^3. To evaluate the skewness with the top-hat smoothing, the integration including 



the window function Wr is tedious and complicated. However, the calculation is tractable in the 



Fourier space because the useful formulae have already been found (Bernardeau 1994a, Bernardeau 
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1994b). Here, utilizing these formulae, we develop the second order perturbation in the Fourier 
space. The variables S m ,5 g and 8 are expanded as follows: 



(x,t) 



d 3 k 



{k,t)e~ ikx , 9(x,t) 



Hereafter, we denote the Fourier coefficient by attaching the hat 



^9(k,t)e- ikx . 



(32) 



The calculation in the Fourier space is thoroughly investigated by Bernardeau, although he 
did not consider the evolution of 5 g ( [Bernardeau 1994~b| ) . We display the results only. The second 
order solutions satisfying the initial conditions are 



4 2) (M) 



d 3 k' 
(2^)3 



1 



Dl(t) -V{k', k-k') + -V{k - k', k') - -C(k', k - k'] 



+-E + (t)C(k',k-k') 



A m (k - k')A m (k'), 



(33) 



5^ (k,t) 5$(k,ti) 
d 3 k' 



+ (D+(t)-l) 



(27T)' 



V(k', k - k')(A m (k')A g (k - k') - A m (k')A m (k - fe')) 



+ 



/; 
6 



(A g (k')A g (k-k>)-(A 2 g )5 D (k>)), 



where 



V(k 1 ,k 2 ) = 1 + 



(ki ■ fe 5 



(ki -k 2 f 



\k\\ 2 ' ' |fci| 2 |A;2| 2 

The solutions (33) and ( |34| ) contain the function E + {t) which satisfies E + (ti 
inhomogeneous solution of the following equation: 

2 + 3 + 



In Einstein-de Sitter universe, we have 



34 

E + (t) = Yl D 2 + (t). 



(34) 

(35) 
1. This is the 

(36) 
(37) 



It is known that the Q and A dependences of E + /D 2 _ is extremely weak ( Bernardeau 1994b| ). 
Therefore, we proceed to analyze the skewness and bi-spectrum by replacing E + (t) with 
(34/21)Z£(t). 



3.1. Skewness 



The second order solutions (33) and ( |34|) together with the linear order solutions give the 
tree-level third order moment by substituting them into (|l6|). From (|i~4|), the result of the 
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smoothed one-point function becomes 
<W> /34 



where 



, , - 7 )(D 2 + -l)(D + -l + b r ) 2 

+ (6 - 7) {b r (D + - 1) + (Z>+ - l) 2 } {(6 r - 1)D + + (b 2 + 1 - 26 r )} 
+ / i -/^{6 r (Z? + -l) + 62} 2 , (38) 



^ [logcr 2 (i?)], (39) 



d(logi?) 

Iw = <?o A f Sl ^ 2 W R {\k x + fc 2 |fl)Ty fl (feiiZ)Ty(fcafl)P(fci)P(fc 2 ). (40) 
The skewness of the galaxy distribution is obtained as follows: 

3,9 = {5 2 {R)) 2 ' tr 2 = - !) + 2 Vo(£>+ - 1) + 6 . (41) 



We see from (41) that S3 j9 evolves from the initial value h ■ Iw and asymptotically approaches the 



skewness of the total mass S3 m , 

34 

S 3 , m = y-7, (42) 

which is constant in time. The numerical value Iw comes from the contribution of non-Gaussian 
initial distribution. It is almost equal to unity for the power spectrum P{k) oc k n (—3 < n < — 1), 
which has been checked by the Monte Carlo integration of (f40|). 

We evaluate the time evolution of skewness by plotting ([II]) for the various parameters. 
Hereafter we assume that the power spectrum simply obeys the single power-law P{k) oc k n , which 
yields 7 = n + 3 from (|39|). In Fig.l, we have plotted the evolution of Ss t9 as a function of the 
expansion factor a. For a specific choice of the parameters, we take = 1, n = —3, h = 3.0, 
for convinience. We have verified that the other choice of the parameters leads to the similar 
behavior of the time evolution of the skewness. The initial conditions at a = 1 for each lines are 
bo = 2.0, ro = 1.0(solid line), bo = 2.0, ro = 0.5(long- dashed line), bo = 1.0, ro = 0.0(short- dashed 
line) and bo = 0.5, ro = 0.0(dotted line). The figure manifests the deviation from the deterministic 
bias r = 1 (solid and long-dashed lines). The skewness with the smaller ro tends to become the 
larger value. We see that the gravitational instability changes the cross correlation largely and 
this enhances the non- linear growth of skewness S^^g. The smaller bias parameter bo also leads 
to the larger skewness (short-dashed and dotted lines). Especially, the skewness in anti-biasing 
case (bo < 1) has a characteristic initial behavior, since the rare distribution of galaxies rapidly 
concentrates on that of total mass due to the attractive force of gravity. 

For more various initial data, the resulting parameters which evolve from the high redshift are 
exhibited in Table. 1. Provided that the initial parameters 60, ro and h(= £3,3) at z = 3, we obtain 
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the bias parameter b and the skewness S^^ g at present, which qualitatively have the same behavior 
as depicted in Fig.l. Although we obtain the skewness Sz g in the treatment of the stochastic bias, 
the non-linear bias parameter b 2 which gives the deterministic 5 g -5 m relation 

S g = f(S m ) = bS m + y (<& - (5 2 J) + ■■■ (43) 

can be evaluated by using the result of the deterministic bias (||). The S g -5 m relation has to do with 
the peak of the joint probability distribution for 5 m and 5 g . The non- vanishing value b 2 represents 
the deviation from a simple linear relation 5 g oc S m . Table. 1 says that the smaller cross correlation 
and the smaller bias parameter bend the 5 g -5 m relation efficiently. The remarkable property of the 
stochastic bias is that 62 could become negative for the small initial cross correlation, although it 
approaches zero asymptotically. This tendency can be interpreted in the following way. Assume 
the existence of the primary linear bias b, it will approach unity asymptotically. However, this 
relaxation may not proceed monotonously. The larger fluctuation accelerates the relaxation 
due to the non-linearity of the gravitational force. Moreover, the small initial cross correlation 
enhances this acceleration. Hence, 5 g -5 m relation with the negative 62 can be obtained. Fry and 
Tegmark & Peebles found that both b and r approach to unity asymptotically (|Fry 1996 , Tegmark 



fc Peebles 1998 ). Here, we found the general tendency of the non-linear bias which is determined 
by the gravitational dynamics. The tilted relation with a negative 62 is consistent with the recent 
numerical simulation ( Dekel and Lahav 1998] ) . 



Table. 1 also gives the cases for the different power spectrum from the scale invariant one 
(n = —3) and the different density parameter from S7o = 1. From the evolution of skewness 
with the spectral index n = — 1, we observe that the initial non- linear bias parameter 62 rapidly 
decreases to zero compared with the case with index n = —3. Since the spectrum with n = — 1 
has larger power than n = —3 on large scales, the result shows that the tilted spectrum could 
enhance the non-linear relation of S g -5 m . As for the low density universe(f2o = 0.2), the growth of 
the density fluctuations becomes slow and the variations of b and b 2 are small. Accordingly, the 
deviations from 6=1 and 62 = in low density universe are larger than that in Einstein-de Sitter 
universe. 



3.2. Bi-spectrum 

In addition to the skewness, the bi-spectrum provides us the other information of three-point 
correlation function. In Fourier space, it is defined by 

(dgi^Sgik^Sgiks)) = 5 D (k X + k 2 + fc 3 )B a (fcl, fe 2 , fc 3 ). (44) 

The amplitude B g relates with the third moments ( |38[ ) as follows: 

(6 3 g (Ro)) = J J j^W(\k 1 + k 2 \R )W(k 1 R )W(k 2 R G )B g (k 1 ,k 2 ,k 3 ). (45) 
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Recall that fci + k 2 + ^3 = 0, we obtain 

4 



{V(k u fe a ) + 7>(fc 2 , fci) - ^(fei, k 2 )}(Dl -!)(£>+- 1 + 6 r ) 2 



B g (ki,k 2 ,k 3 ) 

+ {P(fe!, fc 2 ) + P(fe 2 , fei)}(D+ - 1 + Vo){(Vo - 1)D+ + &g + 1 - 26 r } 
+ ^(& ro(£>+ -l) + ^o) 2 ] P(h)P(k 2 ) 

+ permutations. (46) 



P(/c) is the power spectrum of total mass fluctuations given by (p5[) . A more convenient observable 
is the reduced bi-spectrum amplitude Q g which is briefly discussed in Sec. | (pry 1994| , Fry 1996|) : 



B g (k 1 ,k 2 ,k 3 
P g (h)P g (k 2 ) + P g (k 2 )P g (k 3 ) + P g {kz)P g {k x y 



Qg = TTTTXTTTT-, fTTTTTTTTUT - ( 4 ?) 



The power spectrum of the galaxies P g (k) is written by 

P g (fc) = {( J D+-1) 2 + 2b r (D + -l) + bl) 2 P{k). (48) 

The reduced bi-spectrum amplitude Q m for the mass fluctuation has a remarkable property 
that it is scale independent and constant in time. In our case, Q g is also scale independent due 



to the assumption (|24|), although it depends on time and approaches to Q m . To evaluate the 



configuration shape of the bi-spectrum, we plot Q g as a function of the angle 6 = cos~ 1 (fci • k 2 ) 
by setting k\/k 2 = 2. Fig. 2 shows the various cases of the present amplitude Q g (z = 0) whose 
initial conditions are set at z = 3 in the case of Einstein-de Sitter universe(f2o = 1). We also 
take the initial skewness parameter h = 3.0 and the spectral index n is chosen as n = —3. In 
Fig. 2a, we plot the bi-spectrum for the various initial values of the cross correlation when the 
initial bias parameter is fixed (bo = 2.0). Each solid line shows the result of the initial conditions 
ro = 0.0, 0.25, 0.5, 0.75, and 1.0(/rom up to bottom), respectively. We see that the cross 
correlation affects the shape of the bi-spectrum. For the smaller ro, the shape of Q g becomes 
steeper and it approaches Q m faster. This is the same behavior as the non-linear growth of S*3. 9 . 
In cases the initial cross correlation is fixed, the shape of the bi-spectrum changes significantly 
depending on the bias parameters (Fig. 2b). In the next section, we will show that the small 
shape dependence of the cross correlation plays an important role for explaining both the present 
observation of galaxies and the galaxies at high redshift. 



4. Discussion 

In previous section, we studied the influence of the cross correlation on the time evolution 
of skewness and bi-spectrum. In this section, using these results, we describe how the stochastic 
treatment of the bias affects the analysis of observational data for the galaxy distribution. We will 
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present a method for predicting the bias, cross correlation and skewness at high redshift from a 
present observational data. 

As mentioned in Secffl, the recent observation indicates that the high redshift galaxies at 
z ~ 3 is largely biased, b ~ 6.2 ( peacock 1998| , see also Peacock et al. 1998). On the other hand, 
the correlation functions of the present optical galaxies are one half time larger than that of the 
IRAS galaxies, which means that the bias parameter for the optical galaxies can be evaluated as 



b(z = 0) = 1.5 if we require no bias to the IRAS galaxies ( Hamilton 1997 , Strauss et al. 1992j ). 



This suggests that the time evolution plays an important role for explaining the difference of the 
bias parameter. The present bias parameter is also evaluated by the observation of the higher 
order statistics for various galaxy survey, the Lick catalog, for example. Using the prediction from 
the perturbation theory, Fry examined the deterministic non-linear bias and obtained b = 3.5 and 
&2 = 14.7 to explain the flat shape of Q g for the Lick catalog (Fry 1994). This contradicts with the 
result inferred from IRAS survey. 

In the stochastic description of the bias effect, taking into account the cross correlation 
between 5 m and 8 g , we can fit the bi-spectrum to the Lick catalog so as to be consistent with the 
bias parameter inferred from IRAS galaxy survey. Using the best fit parameters of the stochastic 
bias, we will see the evidence of large biasing for high redshift galaxies and evaluate the initial 
skewness at z = 3. Let us denote 60 an d tq as the initial bias parameter and the initial cross 
correlation set at z = 3, respectively. Providing the present bias parameter as b = 1.5, the initial 
bias parameter 60 is written in terms of the cross correlation tq\ 



^o = 3 Jr 2 Q + 3 - r . (49) 



Substituting the above equation into (|3l|) , b(t) becomes a function of ro- Thus we can know the 
initial bias parameter from the final condition b(z = 0) = 1.5 by varying the ro. Fig. 3 represents 
the time reversal evolution of b for the initial cross correlation < ro < 1. In this inverse problem, 
we have the large bias parameter 60 = 5.20 when ro = 0.0. 

The unknown parameter ro can be determined from the present bi-spectrum of the Lick 
catalog. The initial amplitude and the configuration shape of the bi-spectrum at z = 3 change due 
to the evolution. Using the fact that the spectral index of the Lick catalog becomes n = —1.41 on 
large scales ( Fry 1994| ), the present bi-spectrum Q g (z = 0) can be given as a function of ro and h. 



In Sec.[T^, we observed that the bias parameter 60 affects the configuration shape of Q g (for the 
larger value bo, the shape becomes flatter). In the present case, 60 is associated with ro by (fl9l). 
On the other hand, the initial skewness parameter h can change the amplitude of the bi-spectrum. 
Therefore, adjusting h and ro, we can obtain the bi-spectrum Q g consistent with the Lick catalog. 

In Fig.4, assuming the Einstein-de Sitter universe, the bi-spectrum Q g at z = is depicted by 
varying the initial cross correlation < ro < 1 (solid lines). We use the spectral index n = —1.41 
same as the Lick catalog. The dashed line corresponds to the bi-spectrum of the Lick catalog with 



the best fitting parameters b = 3.5, 62 = 14.7 obtained by Fry (Fry 1994). The figure shows that 
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the configuration shape of the Lick catalog coincides with the bi-spectrum given by the initial cross 
correlation ro = 0.2. Therefore, the bias parameter at high redshift z = 3 is evaluated bo = 4.63 
from Q49|). In this figure, the initial skewness is chosen as h = 6.96. This result indicates that the 
large-scale galaxy distribution significantly deviates from Gaussian distribution and has a large 
bias parameter. This is consistent with the recent observation of Lyman-limit galaxies, which has 
b ~ 6.2 on small scales ~ 7.5Mpc ( peacock 1998[ ). Note that the present skewness evaluated by 
the best fit parameters is £3^ = 4.55 from (|38l ) and (41), which also agrees with the observation 
of the Lick catalog flGaztahaga 1992 ). 



Another explanation of the flat shape of the bi-spectrum in the Lick catalog is that the 
non-linear evolution eliminates the shape dependence ( [Bcoccimarro et al. 1998|) . If this is true, our 
tree-level analysis breaks and the next-to-leading order perturbations (loop correction) are needed 



(Scoccimarro 1997). However, the observation on large scales guarantees the quasi non-linear 
evolution. Thus the most appropriate answer would be the bias effect. We should remark that the 
analysis discussed here can also apply to an arbitrary high redshift galaxies as well as the galaxies 
at z = 3. We can predict the bias parameter, cross correlation and skewness at high redshift if we 
obtain the bias parameter, the cross correlation and the skewness at present consistent with the 
observation of the bi-spectrum. Although this prediction does not take into account the formation 
process of galaxies or morphology of the galaxy type, one can say the presence of cross correlation 
is necessary to lead to the consistent result with the observational data. Therefore, the stochastic 
bias is essential to describe the galaxy distribution. 



5. Conclusion 

In this paper, we have developed the formalism of stochastic bias and investigated the quasi 
non-linear evolution of the galaxy biasing under the influence of gravity. After deriving the 
evolution equations, we analyzed the influence of the cross correlation between 5 m and 5 g on the 
evolution of skewness and bi-spectrum of galaxy distribution. We observed that the stochastic 
description differs from the deterministic biasing. The small value of the cross correlation gives 
rise to the rapid growth of skewness and steep shape of the bi-spectrum. Our main results can be 
summarized as follows: 

• We found the dynamical feature of the stochastic bias from the weakly non-linear analysis. 
The gravity enforces the 5 g -5 m relation to become unity and the larger fluctuation accelerates 
this behavior. The small initial cross correlation enhances the accerelation and leads to 
the non- linear 5 g -5 m relation with a negative curvature 62 < 0. The result implies that the 
gravity significantly affects the stochastic property of the galaxy and the mass. Hence, the 
dynamical evolution of the stochastic bias should be taken into account when we interpret 
the observation of the large scale structure. 

• A method for predicting the bias, cross correlation and skewness at high redshift was 
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presented. Using the present observation of the bi-spectrum and the bias parameter inferred 
from IRAS galaxy survey, the stochastic property of 5 m and 5 g at present can be determined. 
The parameters b, r and £3^ at high redshift are obtainable by the time reversal evolution. 
We solved this inverse problem and predicted the large scale bias b = 4.63 at z = 3 by fitting 
the bi-spectrum to the Lick catalog. The method is easily applicable to an arbitrary high 
redshift survey and will provide a probe for investigating the clustering property of high 
redshift galaxies. 



These results show that the cross correlation plays an important role and the consideration of the 
dynamical evolution of the stochastic bias is required to explain the observational data. Thus, it 
turns out that the stochastic bias and its dynamics are essential to explore the evolution of large 
scale structure. 

Future sky surveys will provide us a lot of information of galaxy distribution. The time 
evolution of the statistical quantities of the large scale distribution will be observed. Since the 
stochasticity of galaxy and mass distribution cannot be removed, to use stochastic bias is inevitable 
for investigating the evolution of large scale structure. Therefore it is important to understand the 
stochastic property of the galaxy and mass distribution by taking into account the time evolution. 
The analytic study of the distribution function may give us the important information for the 
measurement of cosmological parameter and test of inflation theory. As for this direction, we 
have obtained the non- linear 5 m -5 g relation in (8 m , 5 g )-plane although the evaluated relation is a 
kind of expectation value. The most relevant object for the observation is the peaks of the joint 
probability functional which gives the most probable galaxy-mass relation in the measurement of 
galaxy distribution. The primary goal is to clarify the dynamical evolution of the joint probability 
distribution function V[5 m ,5 g ]. Weakly non-linear analysis using the Edgeworth expansion would 
be useful to study the dynamics of the stochastic bias on large scales. Secondary the large scale 
motion of the galaxies should be explored. Since the velocity field of the galaxy distribution 
itself has stochasticity, it is important to examine the effect of stochastic bias on the POTENT 
analysis. Although we studied the relation between 5 m and 5 g in this paper, our formalism should 



be extended to include the velocity field. We will focus on these issues in the future work( Taruya 
Koyama fc Soda 1998; ), 
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Table 1. The present parameters b, b<i and Sz t9 for the various initial conditions given at z = 3 



Initial parameters present values 

(z = 3) (z = 0) 



^0 


n 


bo 


ro 


S 3,g 


b 


b 2 


S 3,g 


1.0 


-3 


0.5 


0.0 


3.0 


0.760 


-0.088 


5.93 


1.0 


-3 


1.0 


0.0 


3.0 


0.790 


-0.133 


5.51 


1.0 


-3 


2.0 


0.0 


3.0 


0.901 


-0.245 


4.48 


1.0 


-3 


2.0 


0.5 


3.0 


1.090 


-0.088 


4.23 


1.0 


-3 


2.0 


1.0 


3.0 


1.250 


0.119 


4.11 


1.0 


-1 


2.0 


0.5 


3.0 


1.090 


0.033 


2.71 


1.0 


-1 


2.0 


0.5 


0.0 


1.090 


-0.128 


2.30 


0.2 


-3 


2.0 


0.5 


3.0 


1.30 


-0.132 


4.49 
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Figure Caption 

Fig.l The time evolutions of the skewness S3 t9 with the specific parameters Qq = 1, n = —3, h = 3.0 
are plotted as a function of expansion factor a. The initial conditions for each lines are 
given at a = 1 as follows: 60 = 2.0, tq = 1.0(solid line), bo = 2.0, tq = 0.5(long-dashed 
line), bo = 1.0, ro = 0.0 (short- dashed line), bo = 0.5, ro = 0.0(dotted line). These lines 
asymptotically approach 34/7, which corresponds to the skewness of the total mass 
distribution. 

Fig. 2 The evolved results of the bi-spectrum amplitude at z = with the spectral index n = —3. 
We plotted Q g as a function of 9 = cos _1 (fci • ^2) setting ki/k% = 2. Each line has the 
same initial skewness h = 3.0 which are given at z = 3 and the universe is assumed to 
be a Einstein-de Sitter universe(Oo = 1): (a). Variation of Q g (6) when the initial bias 
parameter is fixed (60 = 2.0). The initial values of the cross correlation for each solid 
lines are ro = 0.0, 0.25,0.5, 0.75 and 1.0(from top to bottom). We see that the smaller 
cross correlation makes the configuration shape of Q g steeper, (b). Shape and amplitude 
dependences of Q g (9) when the initial cross correlation is fixed (ro = 0.5). The results of the 
bi-spectrum have the initial bias parameters given by 60 = 0.0, 0.5, 1.0, 2.0 and 3.0 (from 
top to bottom). Compared with both figures, influence of the cross correlation is weaker than 
that of the bias parameter. 

Fig. 3 The time reversal evolution of the bias parameters ( fjll) when the present (z = 0) bias 
parameter is given by b = 1.5. The horizontal axis denotes the redshift parameter z. The 
results are shown for the cross correlation ro = 0.0,0.2,0.4,0.6,0.8 and 1.0 (from top to 
bottom), which are the values set at z = 3. The figure shows that the bias parameter takes 
the values from 3 to 5.2 at z = 3. 

Fig. 4 The bi-spectrum Q g by setting the initial conditions at z = 3, all of which have the same 
bias parameter b = 1.5 at present (solid lines). We use the spectral index n = —1.41 same 
as the Lick catalog. Each lines denote the present value of bi-spectrum with the initial cross 
correlation ro = 0.0, 0.2,0.4,0.6,0.8 and 1.0(Jrom top to bottom), respectively. These lines 
have the same initial skewness h = 6.96. The dashed line corresponds to the bi-spectrum of 



the Lick catalog constructed by the time independent non- linear bias ( Fry 1994 ). We can fit 
the bias parameter to the Lick catalog if the initial cross correlation is chosen as ro = 0.2, 
which leads to the initial bias parameter bo = 4.63. This indicates that the large-scale 
distribution of high redshift galaxies is largely biased and has large skewness. 




Fig.l 
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Fig.2a 




Fig.2b 




Fig.3 




Fig.4 



